Transient Effect of Negative Electric Current in Irradiated Semiconductors 
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^P, The peculiarities of electric current are studied occurring in semiconductors with strongly nonuni- 

TlJ ' form distribution of charge carriers. The formation of such nonuniformities and the regulation of 

'^'n ' carrier mobilities can be realized by means of external irradiation, for instance, charge-particle 

, beams and laser irradiation. The transient effect of negative electric current is shown to arise under 

y—^ • some specific conditions. 

^ ! 

> ■ I. INTRODUCTION 

(N ■ 

oo ; 

' The study of electric transport in semiconductors is important for describing and modelling different semiconductor 
00 ■ devices [1,2]. The action of external irradiation can result in the formation in a semiconductor sample of nonuniform 
■ distributions of charge carriers. Thus, ion-beam irradiation leads to the formation of a dense layer of ions located at the 
; distance of their mean free path from the surface. Similar charged layers can be created by irradiating semiconductors 
with other beams of charged particles, say, electrons or positrons. If the injected charges can move, then the irradiated 
material behaves as an extrinsic semiconductor. The mobility of the charge carriers can be activated and regulated by 
involving additionally laser irradiation. Narrow laser beams can also be employed for creating nonuniform distributions 
of charge carriers. 

Transport properties of semiconductors with essentially nonuniform distribution of carriers can be rather specific. 
(~| For instance, in a sample, biased with an external constant voltage, the resulting electric current may turn against 
O . the latter [3,4]. Certainly, this can happen only as a short-time fluctuation after which the current turns back 
O • becoming positive [5,6]. The transient effect of negative electric current has been considered earlier [5,6] for simplified 
^ , models. The aim of the present paper is to give a careful analysis of this effect under conditions typical of realistic 
semiconductor materials. 



II. DRIFT-DIFFUSION EQUATIONS 



Transport properties of semiconductors are usually described by the semiclassical drift-diffusion equations [1,2] 
consisting of the continuity equation 

_+V.j. + -=C., (1) 



with the drift -diffusion current 



and of the Poisson equation 



ji = PiPiE - D^\I , (2) 



eV-^ = 47r(/7i +P2) • (3) 

Here pi = pi{r,t) is a charge density, E — E{f,t) is the electric field, and = ^^(r, t), generation-recombination 
noise [7]. The considered semiconductor is characterized by mobilities p,i, diffusion coefficients Di, and a dielectric 
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permittivity e. In what follows, two types of charge carriers are assumed, positive and negative, such that pi > 
jjLi> Q and p2 < 0, /X2 < 0. The total density of electric current writes 

^ ^ ^ e dE 

3tot=3i+j2 + ^^ , (4) 

being the sum of the drift-diffusion current (2) and the displacement current. 

Let us consider a plane device of area A and width L, so that there is one space variable x G [0, L]. Let the sample 
be biased with an external voltage Vq, which means that 

L 

E{x, t) dx = VQ . (5) 



Define the transit time 



l2 

To=—r-, /i = min{/ii, |/i2|} , (6) 
and also introduce the following characteristic quantities 

Po = , Qo = eAtjQ , tjo = — , 

io=^, D,^„Vo, eo = ^. (7) 
Atq To 

It is convenient to pass to dimensionless notation where the space variable x is measured in units of L; time, in units 
of the transit time (6); and all other quantities in the corresponding units (7). Then the continuity equation (1) is 
written as 

I ^ ( T?\ ^^P^ , Pi c 

^+M.^(p.i?)-A^ + -=C. (8) 

and the Poisson equation (3) becomes 

dE 

— =47r(pi+p2) , (9) 

where the space and time variables are 

< a; < 1 , t > . (10) 

Equation (8) requires an initial and two boundary conditions. The first is given by the initial distribution of charge 
carriers 

Pi(x,0) = /,(x) (z-1,2). (11) 
As the boundary conditions, we may accept the absence of diffusion through the semiconductor surface, which reads 

^p^{x,t) = Q {x = Q,x = l). (12) 

ox 



The role of the boundary condition for Eq. (9) is played by the voltage integral (5) that in the dimensionless notation 

E{x,t)dx=l. (13) 



IS 



/ 

Jo 
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Note that by employing Eqs. (9) and (13), we may write the electric field in the form 



E{x, i) = 1 + 47r 



Q{x,t) — / Q{x,t) dx 
Jo 



(14) 
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where 



Q{x,t)= ( [pi{x' ,t) + p2{x' ,t)\ dx' 
Jo 

As the total current (4), we have 



d \ ( d \ 1 dE 

jtot={f^iE-D^—]p^+ii,2E-D2^]p2 + --^ . (15) 



dxj'^' ' V ^dxj'^^ ' 47r dt 

The quantities whose time behaviour will be of interest for us are the electric current through the semiconductor 
sample 

J{t)= [ jtot{x,t)dx (16) 

and the electric current at the left and right surfaces of the device, 

J{0,t) = jtot{0,t) , J{l,t) = jtot{l,t) . (17) 
The electric current (16), using Eq. (15), can be written as 

J{t)= / [pipi{x,t) + p2P2{x,t)]E{x,t) dx+ 
Jo 

+ £>i[pi(0,i)-pi(l,i)]+L»2[p2(0,t)-p2(l,i)] , (18) 
while for the currents (17) at the left and right surfaces, we get 

J{0,t) = J{t)+ [ [jiQi{x,t)+'y2Q2{x,t)]dx , (19) 
Jo 



where 



and, respectively. 



Q 



r-X ^ 

,{x,t)= pi{x',t)dx', 71 = — (i=l,2), 
^0 n 

J{l,t) = J{0,t)- [ [jipi{x,t)+-f2P2{x,t)]dx . (20) 
Jo 



III. NEGATIVE CURRENT 

In order to demonstrate that there exist conditions when the electric currents defined above can become negative, 
that is, directed against the applied voltage, let us take an illustrative example of narrow initial charge distributions 
which can be approximated by the form 

p,{x,0) = Mx) ^ Q,S{x - ai) . (21) 
Then from Eqs. (18)-(20), we have at the initial time 

J(0) = f^iQiE{ai, 0) + p2Q2E{a2, 0) , (22) 

J(0,0) = J(0)+7igi(l-ai)+72Q2(l-a2) , (23) 

J(l, 0) = J(0) - 7iQiai - 72^202 • (24) 
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For the electric field (14), we get 

E{x, 0) = 1 + AnQi[ai - e(ai - x)] + 47rQ2[ci2 - 6(a2 - x)] , 
where 9(a;) is the unit step function. Using this, for the total current (22), we find 

J(0) = niQi 1 1 + 4nQi (ai - ^) + 47rQ2 K - 6(02 - ai)]| + 

+ l^2Q2 1 1 + 47rQ2 (^a2 - ^ + ^t^Qi [ai - e(ai - 02)] | . (25) 

The relation between the currents (23) and (24) has the form 

J(0,0)- J(l,0)=7iQi+72O2 . 

In what follows, we shall consider two particular cases, when the initial charge distributions are located at the same 
place and when they are separated. 



A. Single— Layer Case 

Assume that both charge distributions (21) are located at the same place 



ai = a2 = a 



The initial electric field, then, is 



E{a, 0) = 1 + AnQ a 



Q = Qi + Q2 



The total current (25) becomes 



J(0) = (M1Q1+M2Q2) 



1 + 47rQ ( a - 2 



Since HiQi > 0, we always have /UiQi + ^2Q2 > 0. Hence, the current (28) can be negative if 

1 1 



a < 



1 

a> - 



2 47rQ 
1 



2 47r|Q| 

As far as < a < 1, inequalities (29) are possible for 

\Q\> 



(Q>o) 
(g < 0) . 



2n 



Inverting Eq. (28), we may define the initial charge location 



1 1 
a = - — 



2 AttQ 



J(0) 



MlQl + A*2(32 



as a function of the current J(0). Measuring the latter gives us the location (31). 
The current (23) at the left surface is negative under the condition 

a (attQ - llRl±J^\ < 27rQ - 1 - ^i^i + ^'^^'^ 
Also, measuring J(0,0), we may define the location 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 
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^ ^ (27rQ - + M2Q2) - (71Q1 + 72Q2) + J(0, 0) .33. 

47rQ(/iiQi + /iaQa) - (71Q1 + 72Q2) 



The current (24) at the right surface becomes negative when 



710^+720^ > 27rQ-l. (34) 



(35) 



MiQi + M2O2 

The location of the initial charge layer can be defined through J(l, 0) as 

(2^Q-l)(AiiQi+M2Q2) + J(1,0) 



A'kQ{hiQi + /U2Q2) - (71Q1 + 72Q2) 



B. Double— Layer Case 



Now consider the case when, at the initial time, two charge layers, described by the distributions (21), are separated 
in space so that 



ai = a < 02 = 1 — a. 
Then, substituting into the electric current (22) the electric fields 

E{ai,Q) = 1 - 27rOi + 4na{Qi - Q2) , 



(36) 



E{a2, 0) = 1 + 2ttQ2 + 47ra(gi - Q2) , 

we have 

J(0) = Mi<3i(l - 27rQi) + M2(32(l + 2t:Q2) + 47ra(Qi - Q2)(AtiOi + At2Q2) • 
This current is negative if 

MiQi - M2Qi 1 



2a{Qi-Q2) < 



Mi<3i + M2Q2 27r 



In particular, when Q2 = —Qi, Eq. (38) yields 



a < 



(37) 



(38) 



4 87rQi 

This inequality, since < a < 5, gives <3i > The conditions for the currents (23) and (24) to be negative are 

71Q1 - 72Q2 



and, respectively, 



MQi - Qi) 



47r(Qi - O2) 



71Q1 - 72Q2 

AtiQi + IJ.2Q2 



MlQl + IJ'2Q2 



AtiQi + IJ.2Q2 



(39) 



(40) 



In this way, we see that each of the electric currents (22)-(24) can become negative at initial time, provided the 
corresponding conditions hold true. 
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IV. QUALITATIVE ANALYSIS 



To understand the general physical picture, we need to solve Eqs. (8) and (9). The generation-recombination noise 
in Eq. (8) is usually modelled by the white Gaussian noise for which the stochastic averaging can be denoted by 
<C . - . This noise is defined by the mean 

<^i(x,t)»=0 (41) 

and by the correlation function 

« ^i{xMj{x',t') » = ^ijS{x - x')5{t - t') , (42) 

where 7y are the parameters characterizing the specific properties of the generation- recombination process. 

An approximate solution of Eqs. (8) and (9) can be found by generalizing the method of scale separation [8-10] to 
the case of differential equations in partial derivatives. To this end, we have to find out which of the functions E, pi 
or p2 could be considered as slow varying and, if so, with respect to what variables. In what follows we assume that 



< pi{x,t) +p2{x,t) > dt 







< oo , (43) 



which will be confirmed a posteriori. Using Eqs. (9) and (43), we get 



1 r d 

lim - / < t)^dt = 0, (44) 



T— >oo T 







hence we can tell that the function E is, on average, slowly varying in space. Then, because of the voltage integral 

(13), we have 



I 



^ d 

<^ —E{x,t) dx = , (45) 
'^t 

which suggests that E is, on average, slowly varying in time. Thus, the function E{x,t) can be treated as a quasi- 
invariant on average, with respect to both space and time. Keeping E fixed in Eq. (8), we obtain a linear equation 
with respect to pi, with constant coefiicients. This equation, complimented by the initial conditions (11) and the 
boundary conditions (12), can be solved. The resulting solution looks a little too cumbersome and we shall not write 
it down here in full, since our aim in this section is to give only a qualitative analysis for understanding the general 
physical picture. Therefore, we shall simplify the solution by considering a thick sample for which, instead of the 
boundary conditions (12), we may formally take 

lim < pj(a;,i) > = . (46) 

X— >±C30 

This simplification is equivalent to passing to an infinite sample, with the simultaneous continuation of pi outside the 
interval [0, 1] by setting pj = for a; < and x > 1. 
Then an approximate solution of Eq. (8) reads 

where the regular part 

/ + 00 
Gi{x-x',t) fi{x')dx' (48) 
-OO 

is caused by the initial distribution fi{x), the Green function being 
and the random part 

pT^'ix, t)= G{x- x', t - t')^i{x', t') dx' dt' (50) 

Jo J-oo 



6 



is generated by the noise £^i{x,t). 

The initial charge distribution can be modelled by the Gaussian 



Zi " 1 2h 

where < aj < 1 and 



/,(a=) = texp|-l^-^^ , (51) 



Qi = y fi{x) dx , J |- ^^ 2b } '^^ ■ 



Strictly speaking, fi{x) is defined as zero for a; < and x > 1. But this restriction can be neglected in the thick-sample 
approximation which, by using the form (51) in the solution (48), yields 

For the random solution (50), because of condition (41), we have 

<p™"(a;,t)>=0. (53) 

And from the definition (42), it follows 



/' 

Jo 



«p™"(ar,%™"(a;',i)» = 

where jij are defined in Eq. (42). From here 

lim < p™"(x, f)p™"(x', f) » = . (55) 

Equations (52) and (53) show that <C ^ exponentially tends to zero as f ^ oo. thus, confirming inequality 
(43). The currents (18)-(20) are influenced by the noise through the correlators (54). The latter, according to Eq. 
(55), are small at short times. Therefore, at the beginning of the process, when t <^ 1, the role of noise is not 
important. This conclusion suggests that for considering transient effects, occurring at f <C 1, the influence of the 
generation-recombination noise may be neglected. 



V. NUMERICAL SOLUTION 



To analyze more accurately the time behaviour of electric current, we have accomplished the numerical calculations 
of Eqs. (8) and (9) with the initial conditions (11) and (51) and the boundary conditions (12) and (13). In agreement 
with the previous section, the noise is neglected. All quantities are given in dimensionless units, as is explained in 
Sec. 2. Varying different parameters entering the problem, we flx /xi = 1 and Qi = 1. We also keep in mind the 
relation D2 = 3Di for the diffusion coefficients, typical of that for holes (Di) and electrons (-D2)- For short, we 
use the notation 71 — 72 = 7 and 61 = 62 = b. Figs. 1 and 2 present the results for the single-layer case, with 
ai = a2 = a; and Figs. 3 and 4, for the double-layer case, when ai = a, 02 = 1 — a. The values of the varying 
parameters are taken so that, when passing to dimensional units, they would correspond to the values characteristic 
for typical semiconductors [1,2]. The current J{t), for Qi > IQ2I, lies always between J(0,t) and J{l,t), so that 
J{l,t) < J{t) < J(0,<), as is clear from Eqs. (18)-(20). Therefore, we concentrate our attention on the behaviour 
of the limiting quantities J(l,t) and J{0,t). The general behaviour of the latter is in agreement with the qualitative 
analysis of Sees. 3 and 4. The principal thing which was impossible to notice in the qualitative analysis is that the 
electric current can become negative not at t = 0, but at some finite time. Anyway, the occurrence of the negative 
current is a transient effect always happening at t <C 1. The dependence of this effect on the physical parameters is 
thoroughly illustrated in the presented figures. 
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Figure Captions 

Fig. 1. Single layer case. 

(a) The electric current J(0, t) at the left surface of a semiconductor sample for the parameters a = 0.25, Di = 
-D2 = 0, 7 = 1, ^2 = —3, for different negative charges: Q2 = (solid line), Q2 = —0.25 (long-dashed line), 
Q2 = —0.5 (short-dashed line), Q2 — —0.75 (dotted line), and Q2 = —1 (dashed-dotted line). 

(a') The electric current J{l,t) at the right surface of semiconductor for the same parameters as in Fig. la. 

(b) The electric current J{0,t) at the left surface for the parameters a = 0.25, 7=1, fJ.2 = —3, Q2 = —0.5, and 
for different diffusion coefScients: £>i = (solid line), £>i = 10""^ (long-dashed line), Di = 10~^ (short-dashed line), 
Di = 10"^ (dotted line). 

(b') The electric current J(l,t) at the right surface for the same parameters as in Fig. lb. 

(c) The electric current J{0,t) at the left surface for the parameters a = 0.25, 7=1, Di = 10'^, Q2 = —0.5, 
and different mobilities of the negative charge carriers: ^2 = —10 (solid line), ^2 = —5 (long-dashed line), ^2 = —3 

(short dashed line). 

(c') The electric current J(l,t) at the right surface for the same parameters as in Fig. Ic. 



Fig. 2. Single-layer case. 

The electric currents at the left, J{0,t) (solid line), and at the right, J{l,t) (long-dashed line), surfaces for the 
diffusion coefficient £>i = 10~^ and different sets of other parameters: 



(a) 


a 


= 0.25, 7 = 


1, m = - 


-3, Q2 = - 


-0.5. 


(b) 


a 


= 0.35, 7 = 


1, = - 


-3, Q2 = - 


-0.5. 


(c) 


a 


= 0.05, 7 = 


10, = 


-3, Q2 = 


-0.5. 


(d) 


a 


= 0.25, 7 = 


10, M2 = 


-3, Q2 = 


-0.5. 


(e) 


a 


= 0.35, 7 = 


10, H2 = 


-3, Q2 = 


-0.5. 


(f) 


a 


= 0.25, 7 = 


0.1, M2 = 


-3, Q2 = 


-0.5. 


(g) 


a 


= 0.25, 7 = 


10, M2 = 


-10, O2 = 


= -0.5 


(h) 


a 


= 0.35, 7 = 


10, M2 = 


-3, Q2 = 


-1. 



Fig. 3. Double-layer case. 
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(a) The electric current J(0, t) at the left surface of a semiconductor sample for the parameters a = 0.1, 7=1, /i2 = 
—3, Q2 = —1, and different diffusion coefficients: Di ^ (solid line), Di = 10^^ (long-dashed line), Di = 10~^ 
(short dashed line), Di = 10~^ (dotted line). 

(a') The electric current J{l,t) at the right surface for the same parameters as in Fig. 3a. 

(b) The left-surface current J(0,i) for the parameters a = 0.1, 7=1, 1^2 = —3, -Di = 10~^, and different initial 
charges: Q2 = (sohd line), Q2 = —0.25 (long-dashed line), Q2 = —0.5 (short-dashed line), Q2 = —0.75 (dotted 
line), and Q2 = —1 (dashed-dotted hne). 

(b') The right-surface current J{l,t) for the same parameters as in Fig. 3b. 

(c) The left-surface current J{0,t) for the parameters a = 0.1, 7=1, Di = 10~^, Q2 = —1, and different 
mobilities: /X2 = —10 (solid line), ^2 = —5 (long-dashed line), /U2 = —3 (short-dashed line). 

(c') The right-surface current J{l,t) for the same parameters as in Fig. 3c. For these parameters the left-and 
right-surface currents are indistinguishable, so the longer time interval is presented here. 

(d) The left-surface current J(0, i) for the parameters a = 0.25, 1x2 = —10, Di = 10^"^, Q2 = —0.1, and different 
relaxation widths: 7 = 25 (solid line), 7 = 10 (long-dashed line), 7=1 (short-dashed line). 

(d') The right-surface current J{l,t) for the same parameters as in Fig. 3d. 



Fig. 4. Double-layer case. 

The electric currents at the left, J(0, t) (solid line), and the right, J(l, t) (long-dashed line), surfaces for the diffusion 
coefficient Di = 10""^ and different sets of other parameters: 

(a) a = 0.1, 7 = 1, = -3, Q2 = -0.5. 

(b) a = 0.1, 7 = 10, = -3, Q2 = -0.5. 

(c) a = 0.1, 7 = 10, H2 = -3, Q2 - -1. 

(d) a = 0.25, 7 = 10, M2 = -3, Q2 = -0.25. 

(e) a = 0.25, 7 = 10, ^2 = -10, Q2 = -0.25. 

(f) a = 0.25, 7 = 10, H2 = -10, O2 = -0.1. 

(g) a = 0.25, 7 = 100, 112 = -10, Q2 = -0.1. 
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